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ABSTRACT 

This  paper  investigates  the  seasonal  variabilities  of  the  South  China  Sea  isopycnal-surface  circulations  and 
of  the  Kuroshio  intrusion  through  the  Luzon  Strait  using  the  U.S.  Navy’s  climatological  temperature  and  salinity 
dataset  (public  domain)  with  Vi  X  V20  resolution  by  the  P- vector  method.  The  representative  pattern  is  a  persistent 
basin-scale  cyclonic  circulation  away  from  the  surface,  and  a  seasonally  varying  circulation  with  a  weak  anti- 
cylonic  gyre  in  the  summer  and  a  strong  cyclonic  gyre  in  the  winter  near  the  surface.  This  pattern  is  consistent 
with  a  classical  view  of  mean  cyclonic  circulation  in  large  stratified  lakes  and  semienclosed  marginal  seas  by 
Emery  and  Csanady  and  with  a  recent  numerical  simulation  using  the  navy’s  Layered  Ocean  Model  by  Metzger 
and  Hurlburt.  The  computed  monthly  volume  transport  through  the  Luzon  Strait  is  negative  (inflow)  all  year 
round  with  a  minimum  value  of  —13.7  Sv  in  February  (strongest  intrusion)  and  a  maximum  value  of  —1.4  Sv 
in  September  (weakest  intrusion).  The  annual  mean  transport  is  —6.5  Sv  (intrusion). 


1.  Introduction 

South  China  Sea  (SCS),  the  largest  marginal  sea  in 
the  West  Pacific  Ocean,  is  separated  from  adjacent 
oceans  by  a  chain  of  islands.  It  contains  a  broad  shallow 
shelf  on  the  south,  extending  along  the  Vietnamese- 
Chinese  coast  to  the  Taiwan  Strait  on  the  north,  a  deep 
sea  basin  in  the  center,  and  a  dangerous  ground  near 
Nansha  with  numerous  reef  islands  over  the  southeast 
(Fig.  1).  It  has  a  bottom  topography  that  makes  it  a 
unique  semienclosed  ocean  basin  overlaid  by  a  pro¬ 
nounced  monsoon  surface  wind.  Extensive  continental 
shelves  (less  than  100  m  deep)  are  found  on  the  western 
and  southern  parts,  while  steep  slopes  with  almost  no 
shelves  are  found  in  the  eastern  part  of  SCS.  The  deepest 
water  is  confined  to  a  bowel-type  trench  and  the  max¬ 
imum  depth  is  approximately  4700  m.  The  SCS  joins 
the  Pacific  Ocean  through  the  Luzon  Strait  between  Tai¬ 
wan  and  the  Philippines.  The  combination  of  geometry, 
connectivity  with  the  Pacific  Ocean,  and  strongly  var¬ 
iable  atmospheric  forcing  contributes  to  the  complex 
dynamics  of  the  flow  in  that  region  (Metzger  and  Hurl¬ 
burt  1996). 

The  Kuroshio,  originating  from  the  North  Equatorial 
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Current,  flows  northward  as  a  western  boundary  current 
east  of  Luzon  and  Taiwan  (Nitani  1970).  The  Luzon 
Strait  is  the  principal  passage  through  which  the  Pacific 
water  enters  the  SCS;  therefore,  the  mass  transport 
through  the  Luzon  Strait  is  critical  in  determining  the 
characteristics  of  SCS  waters  (Shaw  1991).  For  exam¬ 
ple,  Hu  et  al.  (1992)  believed  this  flow  to  contribute  to 
currents  in  the  Taiwan  Strait;  however,  Li  et  al.  (1996) 
showed  that  this  flow  did  not  contribute  to  currents  in 
the  Taiwan  Strait  using  a  numerical  model.  The  volume 
transport  through  the  Luzon  Strait  is  uncertain.  Different 
authors  reported  values  ranging  from  8-10  Sv  by  Huang 
et  al.  (1994)  to  3-Sv  (1  Sv  =  106  m3  s~‘)  outflow  in 
summer  and  3-Sv  inflow  in  winter  by  Wyrtki  (1961). 
Metzger  and  Hurlburt  (1996)  estimated  that  the  mean 
Kuroshio  intrusion  through  the  Luzon  Strait  is  around 
2. 4-4. 4  Sv  (higher  resolution  with  lower  transport  val¬ 
ues)  using  the  U.S.  Navy’s  Layered  Ocean  Model 
(NLOM)  forced  by  wind  stress  climatology  with  various 
resolution. 

The  seasonally  reversing  monsoon  winds  are  also 
play  an  important  role  in  determining  the  upper-ocean 
circulation  (Shaw  and  Chao  1994).  From  April  to  Au¬ 
gust,  the  weaker  southwesterly  summer  monsoon  winds 
result  in  a  wind  stress  of  over  0.1  N  m~2,  which  drives 
a  northward  coastal  jet  off  Vietnam  and  anticyclonic 
SCS  basin-scale  circulation.  From  November  to  March, 
the  stronger  northeasterly  winter  monsoon  winds  cor¬ 
responds  to  a  maximum  wind  stress  of  nearly  0.3  N  m~2 
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Fig.  1.  Geography  and  isobaths  showing  the  bottom  topography  of  the  South  China  Sea. 


causing  a  southward  coastal  jet  and  cyclonic  SCS  basin- 
scale  circulation.  Such  a  seasonal  variability  of  the  sur¬ 
face  circulation  was  simulated  by  several  numerical 
models  (e.g.,  Chu  et  al.  1994,  1996,  1998c,  1999;  Chao 
et  al.  1996).  However,  Metzger  and  Hurlburt  (1996) 
pointed  out  that  the  surface  currents  observed  by  Wyrtki 
(1961)  are  subject  to  Ekman  drift  in  this  highly  variable 
wind  regime. 

Seasonal  occurrence  of  SCS  eddies  have  been  re¬ 
ported  by  several  authors.  Dale  (1956)  and  Uda  and 
Nakao  (1972)  reported  a  cold  eddy  off  the  central  Vi¬ 
etnam  coast  in  summer.  Nitani  (1970)  found  a  cold  eddy 
located  northwest  of  Luzon  in  summer.  Reports  from 
the  South  China  Sea  Institute  of  Oceanology  (1985) 
indicate  that  a  warm-core  eddy  appears  in  summer  and 
winter  in  the  central  SCS,  but  in  summer,  it  is  closer  to 
Vietnam  at  the  surface.  Recently,  a  cold-core  eddy  was 
detected  in  the  central  SCS  during  29  December  1 993— 
5  January  1994,  from  the  analysis  of  TOPEX/Poseidon 
data  (Soong  et  al.  1995).  Chu  et  al.  (1997b)  and  Chu 
and  Chang  (1997)  identified  the  existence  of  a  central 
SCS  surface  warm-core  eddy  in  mid-May  from  a  his¬ 


torical  dataset:  the  U.S.  Navy’s  Master  Observational 
Oceanographic  Data  Set  (MOODS).  From  the  composite 
analysis  on  the  U.S.  National  Centers  for  Environmental 
Prediction  (NCEP)  monthly  sea  surface  temperature 
(SST)  fields  (1982-94),  Chu  et  al.  (1997a)  found  that 
during  the  spring-to-summer  monsoon  transition 
(March-May)  a  warm  anomaly  (greater  than  1.8°C)  is 
formed  in  the  central  SCS  at  1 12°-1 19°30'E,  15°- 
19°30'N.  Emery  and  Csanady  (1973)  reported  that  many 
large  stratified  lakes  and  semienclosed  marginal  seas 
and  estuaries  exhibit  a  persistent  cyclonic  circulation. 
Integrating  NLOM  to  the  statistically  equilibrium  state, 
Metzger  and  Hurlburt  (1996)  simulated  a  persistent  cy¬ 
clonic  circulation  in  the  SCS  upper  layer  with  an  ap¬ 
proximate  thickness  of  250  m. 

Questions  arise:  is  Emery  and  Csanady’s  (1973)  view¬ 
point  or  Metzger  and  Hurlburt’s  (1996)  cyclonic  cir¬ 
culation  pattern  representative  for  the  SCS  circulations 
below  the  mixed  layer?  What  is  the  seasonal  variability 
of  the  Kuroshio  intrusion  through  the  Luzon  Strait?  The 
traditional  geostrophic  calculation  cannot  lead  to  the 
answers  since  only  the  baroclinic  part  is  obtained.  Thus, 
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we  use  the  P-vector  method  (two-step  /3-spiral  method) 
to  invert  the  absolute  geostrophic  circulations  (both  bar- 
otropic  and  baroclinic  parts)  on  the  isopycnal  surfaces 
from  the  navy’s  public  domain  Global  Digital  Environ¬ 
mental  Model  (GDEM)  climatological  monthly  mean 
temperature  and  salinity  data  set  with  V20  X  V20  reso¬ 
lution  (Teague  et  al.  1990). 

The  outline  of  this  paper  is  as  follows.  Section  2 
depicts  the  navy’s  GDEM  dataset.  Section  3  describes 
the  inversion  of  the  absolute  geostrophic  velocity  on  the 
isopycnal  surfaces.  Section  4  depicts  the  thermohaline 
features  at  the  intermediate  depth  using  the  GDEM  data. 
Section  5  describes  the  computation  in  the  isopycnal 
coordinate  system.  Section  6  depicts  the  monthly  mean 
inverted  velocity  fields.  Section  7  discusses  the  seasonal 
variability  of  the  Kuroshio  intrusion  through  the  Luzon 
Strait.  Section  8  presents  the  conclusions. 


(Teague  et  al.  1990).  Different  families  of  representative 
curves  have  been  chosen  for  shallow,  middepth,  and 
deep-depth  ranges,  with  each  chosen  so  that  the  number 
of  parameters  required  to  yield  a  smooth,  mean  profile 
over  the  range  was  minimized.  As  pointed  by  Teague 
et  al.  (1990),  large-scale  oceanographic  features  are  gen¬ 
erally  found  to  be  similarly  represented  in  both  GDEM 
and  the  National  Oceanographic  and  Atmospheric  Ad¬ 
ministration  Climatological  Atlas  of  the  World  Ocean 
temperature  and  salinity.  GDEM  appears  to  render  better 
representations  of  seasonal  variability  and  regions  of 
high  current  shear  because  of  a  different  smoothing 
method  and  a  finer-grid  spacing.  GDEM  data  contains 
the  monthly  mean  temperature  and  salinity  ( T ,  S )  and 
annual  mean  temperature  and  salinity  (T,  S)  fields.  In¬ 
terested  readers  are  referred  to  Teague  et  al.  (1990)  for 
more  informations. 


2.  The  navy’s  GDEM  dataset 

GDEM  is  the  navy’s  global  climatological  monthly 
mean  temperature  and  salinity  data  set  with  a  four-di¬ 
mensional  (latitude,  longitude,  depth,  and  month)  dis¬ 
play.  Data  for  building  the  current  version  of  GDEM 
climatology  for  the  SCS  were  obtained  from  MOODS, 
which  has  144  135  temperature  and  13  768  salinity  pro¬ 
files  for  1930-97.  The  main  limitation  of  the  MOODS 
data  is  its  irregular  distribution  in  time  and  space  (Chu 
et  al.  1997c).  Certain  periods  and  areas  are  over  sampled 
while  others  lack  enough  observations  to  gain  any  mean¬ 
ingful  insights.  Vertical  resolution  and  data  quality  are 
also  highly  variable  depending  much  on  instrument  type 
and  sampling  expertise.  The  monthly  distributions  of 
the  total  SCS  temperature  (Fig.  2a)  and  salinity  (Fig. 
2b)  stations  show  that  the  number  of  temperature  sta¬ 
tions  is  10  times  more  than  the  number  of  salinity  sta¬ 
tions.  The  monthly  distributions  of  the  SCS  temperature 
(Fig.  2c)  and  salinity  (Fig.  2d)  stations  with  the  profile 
deeper  than  400-m  depth  are  approximately  one  third 
of  the  total  stations.  Although  the  number  of  observa¬ 
tions  below  400  m  is  much  smaller,  the  monthly  station 
distributions  (Figs.  2c, d)  still  show  sufficient  data  to 
build  up  the  GDEM  climatology. 

Yearly  temperature  (Fig.  3a)  and  salinity  (Fig.  3b) 
profile  numbers  show  temporally  uneven  distribution 
with  almost  no  observations  in  the  whole  SCS  in  certain 
years  (e.g.,  1944  for  temperature,  and  1944-46,  1952— 
54,  1993-95  for  salinity)  and  many  observations  in  other 
years  (e.g.,  more  than  12  000  temperature  profiles  in 
1966  and  1968,  and  more  than  1100  salinity  profiles  in 
1981).  Spatial  and  temporal  irregularities  along  with  the 
lack  of  data  in  certain  regions  must  be  carefully  weight¬ 
ed  in  order  to  avoid  statistically  induced  variability. 

The  basic  design  concept  of  GDEM  is  the  determi¬ 
nation  of  a  set  of  analytical  curves  that  represent  the 
mean  vertical  distributions  of  temperature  and  salinity 
for  grid  cells  (0.5°  X  0.5°)  through  the  averaging  of  the 
coefficients  for  the  curves  found  for  individual  profiles 


3.  Determination  of  absolute  geostrophic  flow 

Chu  et  al.  (1998c)  showed  from  a  numerical  study 
that  the  oceanographic  conditions  of  SCS  is  closely  re¬ 
lated  to  monsoon,  and  especially  the  wind  stress  largely 
affect  the  near-surface  (30-m  depth)  seasonal  circula¬ 
tion.  However,  the  wind  effect  weakens  with  the  depth. 
The  thermohaline  circulation  becomes  important  in  the 
sub-surface  layer.  In  this  study,  we  use  the  GDEM  data 
to  invert  the  absolute  geostrophic  velocity  on  isopycnal 
surfaces.  The  inverse  method  proposed  in  this  section 
is  the  two-step  /3-spiral  method:  (a)  determining  the  di¬ 
rection  of  the  velocity,  and  then  (b)  determining  the 
magnitude  of  the  velocity. 


a.  Geostrophic  flow  on  isopycnal  coordinates 

As  pointed  out  by  Stommel  and  Schott  (1977)  and 
Wunsch  (1978),  in  determining  large-scale  circulation 
from  hydrographic  data,  we  can  be  reasonably  confident 
on  the  assumptions  of  geostrophic  balance,  mass  con¬ 
servation,  and  no  major  cross-isopycnal  mixing  (except 
water  masses  are  in  contact  with  the  atmosphere).  The 
potential  density  (p)  of  each  fluid  element  would  be 
conserved. 

The  isopycnal  coordinate  system  (x,  y,  p)  is  used  with 
x  axis  in  the  zonal  direction,  y  axis  in  the  longitudinal 
direction,  p  axis  in  vertical.  The  geostrophic  balanced 
motion  on  an  isopycnal  surface,  p,  with  a  pressure,  p, 
is  given  by  (Bleck  and  Smith  1990) 


fu 


m  m 

— >  fv  =  — , 

dy  dx 


(1) 


where  u,  v  are  velocity  components  in  x  and  y  axes,  M 
=  pip  +  gz,  is  the  Montgomery  potential.  The  hydro¬ 
static  balance  is  written  by 


p  9M 
—  +  —  =  0. 
p-  dp 


(2) 
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Fig.  2.  Spatial  distributions  of  MOODS  stations  during  1930-97:  (a)  total  temperature  stations,  (b)  total  salinity  stations,  (c)  temperature 
stations  with  depths  deeper  than  400  m,  and  (d)  salinity  stations  with  depths  deeper  than  400  m. 
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Fig.  2.  ( Continued ) 
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Fig.  2.  ( Continued ) 
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Fig.  3.  Temporal  distributions  of  MOODS  stations  during  1930-97:  (a)  temperature  and  (b) 

salinity. 


The  adiabatic  density  conservation  and  incompressible  b.  Tangential  unit  vector 
continuity  equations  are  combined  to  give  an  equation 

for  dp/dp.  For  the  steady  state,  the  equation  is  From  (3)  we  define  a  streamfunction  4'(x,  y)  by 


Afe,)  +  A(AO  =  o. 

dx\dp  J  dy\dp  ) 


(3) 


Notice  that  the  differentiations  with  respect  to  x  and  y 
are  on  the  isopycnal  surface.  Equations  (l)-(3)  are  the 
basic  equations  for  the  inversion. 

Differentiation  of  (1)  with  respect  to  p  and  use  of  (2) 
lead  to  the  thermal  wind  relation 


q  d'F  q  dV 

- ,  fv= - , 

g  dy  g  dx 


where  g  is  the  gravitational  acceleration  and 


q  =  fg 


dp\ 

dp) 


-1 


(5) 


(6) 


dll  l  dp  dv  1  dp 

dp  fp2  dy’  dp  fp2dx' 


is  the  potential  vorticity  in  the  isopycnal  coordinate  sys- 
(4)  tem.  Neglect  of  relative  vorticity  in  the  q  computation 
may  induce  a  small  but  systematic  error  into  estimation 
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Fig.  4.  Isoline  of  potential  vorticity  ( q )  is  the  trajectory  on  the  iso- 
pycnal  surface. 


of  potential  vorticity  (G.  T.  Needier  1986,  unpublished 
manuscript).  From  (1)  and  (5)  we  have 

d'jr/dx  _  dW/dy 
dM/dx  ~  dM/dy  ’ 

which  implies  that  for  a  given  p  the  streamfunction  4' 
is  a  function  of  the  Montgomery  potential  M  only,  4' 
=  4TM).  Thus,  Eq.  (5)  becomes 


/« 


g 


m 

dy  ’ 


fy 


Comparing  (7)  to  (1)  we  have 


4' TV/) 

g 


9M 

dx 


(7) 


q  = 


g 

4 r'(M) 


=  q(M) 


which  implies  that  for  a  given  p  the  potential  vorticity 
q  is  a  function  of  M  only,  and  in  turn  leads  to  the 
conservation  of  potential  vorticity  on  the  isopycnal  sur¬ 
face. 


dq  dq 

u - r  v —  =  (J. 

dx  dy 


(8) 


The  potential  vorticity  conservation  leads  to  the  fact 
that  any  water  particle  should  move  along  q  isolates  on 
the  isopycnal  surface,  that  is,  any  q  isoline  is  a  trajectory 
of  water  particles  (Fig.  4).  For  each  trajectory,  we  can 
define  a  unit  tangential  vector. 


_  dq\ 
IV^lydy  dx^ ) 


(9) 


Such  a  two-dimensional  vector,  c,  is  equivalent  to  the 
earlier  defined  P  vector  on  z  levels  (Chu  1995;  Chu  et 
al.  1998a,b).  The  necessary  condition  for  the  existence 
of  the  vector  c  is 


Vq  ^  0.  (10) 

At  any  point  on  the  isopycnal  surface,  the  vector  c 
=  ( cx ,  cy)  indicates  the  tangential  direction  of  the  tra¬ 
jectory  and  therefore,  is  parallel  to  the  absolute  velocity 
vector, 

V  =  y(x,  y,  p)c,  (11) 

where  the  absolute  value  of  y  is  the  speed,  lyl  =  I VI. 


c.  Determination  of  the  proportionality  y 

Applying  the  thermal  wind  relation  (4)  to  any  two 
different  isopycnal  levels  pk  and  pm,  a  set  of  algebraic 
equations  for  determining  the  parameter  r  is  obtained: 

y(<r)c(«  _  yWcW  =  Am ^ 

y®cf  -  =  Avkm,  (12) 

which  are  two  linear  algebraic  equations  for  y{k)  and 
y(m).  Here  yu)  =  y(x,  y,  p,),  and 


(Am*,,,,  AvJ 


I  [Pk  1 JdJL 

f  L  p2w 


dp. 


As  soon  as  y®  is  obtained,  the  velocity  field  can  be 
computed  by  (11). 

If  the  determinant 


'-x 

c»> 


^  0, 


(13) 


the  algebraic  equations  (12)  have  definite  solutions  for 
ym  (m  k)\ 


Am^ 

r(m) 

^ X 

g 

<1 

r(m) 

.(*) 

X 

n  (w) 

.(*) 

V 

r  (m) 

(14) 


Now,  the  velocity  inversion  can  be  divided  into  two 
steps:  (i)  determination  of  the  unit  vector  c,  and  (ii) 
determination  of  y.  Chu  et  al.  ( 1998a, b)  show  the  ad¬ 
vantage  of  this  treatment  in  the  velocity  inversion.  For 
water  columns  satisfying  the  two  necessary  conditions 
(10)  and  (13),  we  may  use  (14)  to  compute  y{k>  for  the 
level  pk.  There  are  N  —  1  sets  (m  =  1,2 ,  k  —  1,£  + 

l,  .  .  . ,  N)  of  equations  (12)  for  calculating  y(k).  Here 
N  is  the  total  number  of  vertical  levels  of  the  water 
column.  All  the  N  —  1  sets  of  equations  are  compatible 
under  the  thermal  wind  constraint  and  should  provide 
the  same  solution.  However,  because  of  errors  in  mea¬ 
surements  (instrumentation  errors)  and  computations 
(truncation  errors),  the  parameters  y®  may  vary  with 

m.  A  least  squares  error  algorithm  was  used  to  minimize 
the  error.  For  further  details  and  validation  of  the  al¬ 
gorithm  presented  here  see  Chu  et  al.  (1998a, b).  (In¬ 
terested  readers  can  obtain  the  software  by  contacting 
Mr.  C.  Fan  at  fan@nps.navy.mil  or  online  at  http:// 
web.nps.navy.mil/~chu.) 


4.  Monthly  mean  thermohaline  fields 

The  monthly  mean  three-dimensional  temperature 
and  salinity  fields  obtained  from  the  GDEM  dataset  is 
similar  to  the  climatological  monthly  mean  fields  di¬ 
rectly  computed  from  the  MOODS  as  depicted  in  Chu 
et  al.  (1999).  The  sea  surface  temperature  (SST)  vari- 
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ation  obtained  from  the  GDEM  data  is  quite  consistent 
with  earlier  investigations  based  on  NCEP  data  (Chu  et 
al.  1997a)  and  based  on  the  MOODS  data  (Chu  et  al. 
1997b,  1999).  It  is  not  our  intention  to  present  a  detail 
SCS  climatological  thermohaline  structure  in  this  paper. 
We  only  show  some  representative  features  at  150-m 
depth  as  illustration. 

a.  Temperature 

Figure  5  is  the  horizontal  depiction  of  monthly  mean 
temperature  at  150-m  depth.  The  contour  interval  is 
0.5°C.  In  them,  we  see  that  the  SCS  basin  is  dominated 
by  a  cool  pool  with  several  cold  centers  located  north¬ 
west  of  Luzon  and  south  of  Dongsha  (DS),  near  South 
Vietnamese  Bight  (SVB),  Xisha  (XS),  and  Liyue  Bank. 
All  of  these  cold  centers  were  previously  identified  by 
an  airborne  expandable  bathythermograph  survey  (Chu 
et  al.  1998a)  and  statistical  analysis  on  the  MOODS 
data  (Chu  et  al.  1997a,  1999).  However,  the  individual 
cold  center  was  observed  much  earlier.  For  example. 
Dale  (1956)  and  Uda  and  Nakao  (1972)  reported  a  cold 
eddy  off  the  central  Vietnam  Bight  in  the  summer  (i.e., 
the  SVB  cold  center).  Nitani  (1970)  found  a  cold  eddy 
located  northwest  of  Luzon  in  summer  (DS  cold  center). 
The  summer  SVB  and  DS  cold  centers  were  also  sim¬ 
ulated  by  Chao  et  al.  (1996)  using  a  primitive  equation 
model.  The  warm  water  appears  in  the  vicinity  of  the 
SCS  cool  water.  A  transient  central  SCS  warm  pool  was 
identified  in  spring  (Chu  et  al.  1997a, b).  A  warm  water 
mass  is  situated  in  the  Philippine  Sea  all  year  round  east 
of  the  Luzon  Strait,  where  a  strong  front  appears  be¬ 
tween  the  warm  Kuroshio  water  and  cool  SCS  water. 
Other  two  warm  centers  appear  southwest  of  Luzon 
(SWL),  and  northwest  of  Borneo. 

The  SCS  seasonal  thermal  variability  can  be  char¬ 
acterized  by  strengthening/weakening  of  the  major  SCS 
cool  pool.  The  basin-wide  cool  pool  expands  during  the 
winter  and  shrinks  during  the  summer.  In  December, 
this  cool  pool  occupies  the  vast  area  of  SCS  from  9°  to 
20°N  with  two  evident  cold  centers:  the  Northwest  Lu¬ 
zon  cold  center  (located  at  116°-120°E,  16°-19°N) 
strengthens  during  the  winter  monsoon  season  and  how¬ 
ever,  the  Xisha  cold  center  weakens  and  disappears  in 
February  and  March.  On  the  other  hand,  the  SVB  cold 
center  (located  at  109°-112°E,  12°-16°N)  appears  in 
May  and  strengthens  during  the  summer  monsoon  sea¬ 
son,  and  at  the  same  time  the  Northwest  Luzon  cold 
center  weakens  and  disappears. 

The  warm  tongue  penetrates  northwestward  right  af¬ 
ter  passing  through  the  Luzon  Strait  from  the  Kuroshio 
and  then  stretches  southwestward  along  the  western 
boundary  of  the  SCS  with  a  very  narrow  width  (100- 
200  km).  This  is  consistent  with  earlier  research  on 
climatological  data  (Shaw  1989,  1991)  and  in  situ  ob¬ 
servational  data  (Li  et  al.  1998). 

The  SWL  warm  center  (located  at  116°-120°E,  12°- 
14°N)  occurs  in  March.  The  warm  tongue  associated 


with  the  SWL  warm  center  expands  westward  during 
the  spring-to-summer  transition,  and  occupies  the  east 
part  of  central  SCS  basin  in  the  summer  monsoon  season 
(June-September).  During  the  autumn-to-winter  tran¬ 
sition  (October-November)  the  XS  cold  center  occurs 
and  strengthens  in  December. 

b.  Salinity 

Figure  6  is  the  horizontal  depiction  of  monthly  mean 
salinity  at  150-m  depth.  The  contour  interval  is  0.05 
psu.  The  most  evident  feature  is  the  advance/retreat  of 
the  salty  tongue  represented  by  a  34.60-psu  isoline  from 
the  Philippine  Sea  through  Luzon  Strait.  This  salty 
tongue  penetrates  northwestward  right  after  passing 
through  the  Luzon  Strait  and  then  stretches  southwest- 
ward  along  the  western  boundary  of  the  SCS.  The  salty 
tongue  has  the  minimum  penetration  into  the  SCS  in 
October.  As  the  salt  tongue  advances  into  the  north¬ 
western  SCS  along  the  Chinese-Vietnamese  coast,  a 
salty  center  occurs  near  Dongsha  (1150-117°E,  19°- 
21°N)  in  December.  As  the  salty  tongue  further  advances 
in  the  next  six  months  (January-June),  the  Dongsha 
salty  center  weakens  and  disappears  in  March,  and  an 
evident  salty  center  shows  up  near  the  SVB  (110°- 
113°E,  14°-16°N)  for  six  months  (January-June)  with 
salinity  greater  than  34.7  psu  in  February,  March,  April, 
and  greater  than  34.65  psu  in  other  months.  June  is 
characterized  as  the  maximum  penetration  of  the  salty 
tongue  into  the  SCS.  The  salty  tongue  starts  to  retreat 
in  July  and  the  SVB  salty  center  weakens.  A  salty  center 
appears  in  southeastern  SCS  west  of  Palawan-Borneo 
(1 14°-1 16°E,  8°-ll°N)  enclosed  by  34.6  psu.  This  salty 
center  expands  toward  the  west  until  December. 

Such  a  feature  is  somewhat  inconsistent  with  the  ear¬ 
lier  studies  (e.g.,  Shaw  1991).  This  is  due  to  the  sparse 
salinity  observations  in  the  SCS.  The  MOODS  dataset 
for  building  the  GDEM  climatology  has  much  less  sa¬ 
linity  than  temperature  stations  (Fig.  2).  The  seasonal 
variability  depicted  here  is  not  reliable  for  certain  re¬ 
gions  such  as  the  region  between  Liyue  bank  and  Nan- 
sha,  where  there  is  no  observational  data. 

5.  Computation  in  the  isopycnal  coordinate  system 

a.  Vertical  discretization 

Following  earlier  studies  (e.g.,  Keffer  1985;  Tally 
1988),  we  use  three  reference  levels  for  potential  density 
computation:  cre  for  the  surface,  er2  for  2000  decibar 
(db),  and  cr4  for  4000  db.  Use  of  <x2  and  cr4  provides  a 
better  potential  vorticity  representation  for  levels  near 
2000  and  4000  db,  respectively.  Here, 

irm  =  pm  —  1000  kg  m~3,  m  =  0,  2,  4,  (15) 

where  cr0  =  cre.  The  potential  density  cre,  cr2,  and  cr4 
represent  upper,  intermediate,  and  deep  layers,  respec¬ 
tively.  We  consider  the  following  ranges  for  the 
crm-values: 
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Fig.  5.  Monthly  mean  temperature  (°C)  fields  at  150-m  depth  obtained  from  the  GDEM. 
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Fig.  6.  Monthly  mean  salinity  (psu)  fields  at  150-rn  depth  obtained  from  the  GDEM. 
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22.200  <  ag  <  27.725,  31.625  <  a2  <  37.150, 

41.30  <  tr4  <  45.90.  (16) 

Notice  that  the  values  of  27.725,  37.15,  45.90  are  max¬ 
imum  values  for  cr0,  cr2,  and  cr4  computed  from  the 
GDEM  T,  S  dataset.  We  further  discretize  cr0,  er2,  and 
cr4  with  the  same  increment, 

A  a  =  0.025  kg  m-3,  (17) 

for  cr0  and  cr2,  and 

Act  =  0.02  kg  m~3 

for  cr4.  Thus,  we  have  222  rr„- layers,  222  cr2-layers,  and 
280  <x4-layers.  Within  each  layer,  the  density  is  vertically 
uniform. 

In  order  to  well  resolve  isopycnal  surfaces,  we  use 
cubic  spline  to  interpolate  the  T,  S  data  into  246  z  levels 
with  three  different  increments:  5  m  from  0-  to  100-m 
depth,  10  m  from  100-  to  1000-m  depths,  20  m  from 
1000-  to  2500-m  depths,  and  50  m  below  2500-m 
depths.  Thus,  we  build  up  a  high  resolution  /-coordinate 
data  set  [T{zf),  S(Zj),  <x„(zJ)].  The  symbol  overcarat  in¬ 
dicates  either  the  observational  data  or  the  values  com¬ 
puted  from  the  observational  data.  Notice  that  the  values 
of  27.725  is  the  maximum  value  for  cre,  computed  from 
the  GDEM  T,  S  dataset.  For  simplicity  and  no  loss  of 
generality,  we  will  use  er„  for  illustration. 


b.  Transformation  of  T,  S  data  from  a  z  to  cre 
coordinate  system 

The  transformation  is  fulfilled  by  comparing  the 
z-coordinate  potential  density  data  &g{Zj)  with  the  dis¬ 
crete  cr8-values  at  the  bottom  of  the  kth-isopycnal  layer, 
o-bf,k), 

abg{k)  =  22.1875  +  (k  -  l)Acre, 
c rg(k )  =  ^K(k)  +  ag(k  +  1)].  (18) 


The  geometric  depth  for  the  bottom  of  the  cre(k)-layer 
is  obtained  by 


Dbk  = 


Dbk  =  -Zj  ~ 


zp  if  &g(Zj)  =  crg(k)  and 

-  <T,(z,) 


(19) 


(Zj+i  -  Zj), 


&e(Zj+ 1)  -  a-„(Zj)"’J+' 

if  o-s(Zj)  <  o-b(k )  <  &g(zj+1),  (20) 


where  the  superscript  b  indicates  the  bottom  of  the  kth- 
isopycnal  layer.  The  thickness  of  the  kth-isopycnal  layer 
is  obtained  by 

K  =  Db-  DU.  (21) 


c.  Pressure  and  potential  vorticity  computation 

After  hk  is  obtained,  we  may  compute  the  pressure 
and  potential  vorticity  fields.  Starting  from  the  surface 


I 

I 

I 


(k  =  1),  the  kth  layer  has  the  density  pk,  and  thickness 
hk,  as  shown  in  Fig.  7.  Since  the  pressure  is  not  uniform 
within  the  layer,  it  has  the  value  of 

P'k  =  g  'n  pfii  (22) 

1=1 

at  the  top  of  the  kth  layer,  and  the  value  of 

Pbk  =  gli  pfii  (23) 

1=1 

at  the  bottom  of  the  kth  layer.  We  may  use  the  average 


pi  +  pb  (U  ~  1  „  \ 

Pk  =  — ^ -  =  2  Pihi  +  2Pkhk)  (24) 

to  indicate  the  pressure  at  the  middle  of  the  kth  layer. 
For  the  kth-isopycnal  layer  cre(k),  the  potential  vorticity 
(6)  is  discretized  by 


qk  = 


fgS<Te 

Pk  ~  Pk 


(25) 


Earlier  work  (e.g.,  McCartney  1982;  Keffer  1985;  Talley 
1988)  also  shows  the  benefit  of  using  potential  vorticity 
in  ocean  circulation  studies. 


6.  Monthly  mean  circulations 

Most  important  features  from  our  computation  are  the 
Kuroshio  intrusion,  basin  gyre,  flow  separation,  and  me- 
soscale  eddies.  The  inverted  monthly  SCS  circulation 
agrees  quite  well  with  earlier  observational  study  (Wyrt- 
ki  1961)  at  the  surface  and  with  modeling  results  (Chu 
et  al.  1999)  at  the  deeper  levels. 
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a.  Subsurface  (crg  =  23.0  kg  m  3) 

The  subsurface  shows  a  strong  seasonal  variation 
(Fig.  8).  Due  to  the  shallowness  of  <xe  =  23,  the  Ku- 
roshio  east  of  the  Luzon  Strait  is  not  well  inverted. 
Warm  and  salty  Kuroshio  water  enters  the  SCS  through 
the  Luzon  Strait  evidently  in  October-March,  the  tran¬ 
sition  to  and  during  the  winter  monsoon  season.  The 
winter  circulation  pattern  is  the  basin-wide  cyclonic 
gyre.  A  cyclonic  mesoscale  eddy  splits  from  the  gyre 
near  DS  in  February,  strengthens  afterward,  and  be¬ 
comes  evident  through  the  whole  summer  monsoon  sea¬ 
son.  The  Kuroshio  intruded  flow  is  likely  to  be  confined 
near  the  continental  slope  south  of  China.  The  DS  cy¬ 
clonic  eddy  is  undetectable  from  November  to  January. 
On  the  other  hand,  the  basin-wide  cyclonic  gyre  weak¬ 
ens  in  April  with  the  disappearance  of  the  southward 
western  boundary  currents  near  the  Vietnamese  coast. 

The  summer  (June-August)  pattern  is  generally  a 
weak  anticyclonic  gyre.  A  northward  western  boundary 
current  occurs  from  the  south  to  the  Vietnamese  coast 
with  a  much  weaker  speed  (less  than  0.1  m  s_l)  than 
its  counterpart  (southward  western  boundary  current)  in 
the  winter.  The  western  boundary  current  splits  into  two 
currents  at  10°N:  the  coastal  current  and  the  offshore 
current.  The  offshore  current  recirculates  and  forms  a 
mesoscale  anticyclonic  eddy  near  Natuna  Island  (NI). 
The  NI  eddy  was  simulated  using  a  numerical  model 
(Chu  et  al.  1999).  The  coastal  branch  continues  north 
then  east  at  Hainan  Island. 

b.  Intermediate  level  (crB  =  26.2  kg  m  3) 

As  the  water  depth  increases,  the  seasonal  variation 
of  the  inverted  velocity  reduces  with  depth.  The  inter¬ 
mediate  level  (cre  =  26.2),  changing  from  280  to  460 
db,  has  a  persistent  basin-scale  cyclonic  gyre  as  the 
major  feature  of  the  circulation.  There  is  no  seasonal 
reversal  of  the  western  boundary  current.  The  seasonal 
circulation  variability  is  features  by  strengthening, 
weakening,  and  splitting  of  a  basin-wide  cyclonic  gyre 
(Fig.  9).  This  gyre  is  evident  and  strong  in  the  winter 
and  weakens  in  the  summer.  As  spring  approaches 
(April-May),  the  western  boundary  currents  (south¬ 
ward)  near  the  SVB  weakens.  This  “broken”  gyre  fea¬ 
ture  continues  in  the  whole  summer  (June-August).  The 
western  boundary  currents  (southward)  intensifies  in  the 
fall.  The  cyclonic  basin-scale  gyre  all  year  round  at  the 
intermediate  level  consists  with  the  numerical  models 
(Metzger  and  Hurlburt  1996;  Chu  et  al.  1999). 

7.  Kuroshio  intrusion  through  the  Luzon  Strait 

The  seasonal  variation  of  the  intrusion  of  the  Kuro¬ 
shio  Water  into  the  SCS  through  the  Luzon  Strait  has 
been  investigated  in  earlier  studies  (Shaw  1989,  1991). 
As  pointed  out  by  Shaw  (1991),  the  intrusion  process 
is  important  not  only  to  the  transport  of  water  masses 


into  the  SCS,  but  also  to  the  shelf  circulation  off  the 
southeast  coast  of  China.  Shaw  (1989)  used  the  dis¬ 
criminant  analysis  method  to  classify  the  water  mass  T, 
S  characteristics  at  150,  200,  and  250  m,  and  found  that 
water  characteristics  of  the  Philippine  Sea  (Kuroshio) 
were  identifiable  along  the  continental  margin  south  of 
China  from  October  to  January.  The  presence  of  this 
water  indicated  an  intrusion  current  from  the  Philippine 
Sea  into  the  SCS. 

We  use  the  east-west  geostrophic  velocities  on  iso- 
pycnal  surfaces  along  120°45'E  longitude  with  the  min¬ 
imum  distance  between  Taiwan  and  Luzon  to  represent 
the  flow  across  the  Luzon  Strait  (Fig.  10).  Positive  val¬ 
ues  show  the  outflow  from  the  SCS  to  the  Philippine 
Sea,  and  negative  values  represent  the  inflow  from  the 
Philippine  Sea  to  the  SCS  (Kuroshio  intrusion).  The 
estimate  of  the  transport  through  the  Luzon  Strait  is 
sensitive  to  the  cross-section  selection  due  to  larger  gra¬ 
dient  of  the  property  (temperature,  salinity,  and  Mont¬ 
gomery  potential,  etc.)  in  the  zonal  direction  than  in  the 
latitudinal  direction.  The  signal  of  the  zonal  flow  might 
be  contaminated  by  the  fluctuation  of  the  Kuroshio. 

Flow  through  the  Luzon  Strait  has  a  strong  seasonal 
variation  in  the  upper  and  intermediate  layers.  Evident 
Kuroshio  intrusion  occurs  from  November  to  April  from 
the  surface  to  cre  =  27.0  with  the  maximum  speed  of 
0.18  m  s_1  in  the  middle  of  the  Strait  at  cr„  =  23.75  in 
December.  From  January  to  March  the  core  of  this  jet 
splits  into  two  jet  cores.  The  upper  jet  core  (cr8  =  24) 
shifts  toward  the  south  (Luzon  Island)  with  the  maxi¬ 
mum  speed  reducing  from  0.17  m  s~'  (January)  to  0.14 
m  s_1  (March).  The  lower  jet  core  (cre  =  25)  keeps  the 
same  location  in  the  middle  of  the  strait  (20°-20°30'N). 
In  the  spring  and  summer,  the  Kuroshio  intrusion  weak¬ 
ens  drastically.  The  maximum  speed  becomes  0.06  m 
s  1  at  the  intermediate  level  (cre  =  26.0)  in  the  southern 
part  of  the  Luzon  Strait  (near  19.5°N)  from  May  to 
October.  Weak  outflow  (0.02  m  s  1 )  to  the  Philippine 
Sea  occurs  in  the  northern  part  of  the  Luzon  Strait  (north 
of  20°45'N)  from  the  surface  to  cre  =  25.5  during  the 
summer  monsoon  season. 

We  computed  monthly  total  and  layered  (between  two 
ue  levels)  volume  transports  through  the  Luzon  Strait 
(Fig.  11).  The  dashed,  dotted,  and  solid  curves  indicate 
outflow  (eastward),  inflow  (westward),  and  net  trans¬ 
ports.  The  outflow  is  much  weaker  than  the  inflow  for 
the  total  transport  and  the  most  layered  transports  except 
the  deep  layer  (cre  >  27.025),  where  the  outflow  is  dom¬ 
inant.  The  total  transport  is  negative  (inflow)  all  year 
round  with  a  minimum  value  of  —13.7  Sv  in  February 
and  a  maximum  value  of  —1.4  Sv  in  September  with 
the  annual  transport  of  —6.5  Sv;  which  is  between  Metz¬ 
ger  and  Hurlburt’s  (1996)  estimation  (2. 4-4.4  Sv)  and 
Huang  et  al.’s  (1994)  estimation  (8-10  Sv).  The  total 
and  layered  volume  transports  at  the  Luzon  Strait  might 
be  useful  for  SCS  numerical  modeling. 

There  is  a  time  lag  in  seasonal  variations  between  the 
Kuroshio  intrusion  and  the  salty  tongue  penetration  (Fig. 
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Fig.  8.  Monthly  mean  subsurface  (<r9  =  23.0  kg  m  3)  velocity  vector  field. 
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Fig.  9.  Monthly  mean  intermediate  level  (<r8  =  26.2  kg  m  3)  velocity  vector  field. 
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Fig.  10.  Monthly  mean  east-west  geostrophic  velocities  (eastward  positive)  on  isopycnal  surfaces  along  120°45'E  longitude  representing 

the  flow  across  the  Luzon  Strait  (unit:  0.01  m  s_1). 


2436 


JOURNAL  OF  PHYSICAL  OCEANOGRAPHY 


Volume  30 


Fig.  1 1 .  Monthly  total  and  layered  (between  two  ag  levels)  volume  transports  through  the  Luzon  Strait 
along  120°45'E  longitude:  (a)  total,  (b)  cre  less  than  24.0,  (c)  crg  from  24.0  to  25.0,  (d)  crg  from  25.0  to  26.0, 
(e)  ag  from  26.0  to  27.0,  and  (f)  ag  more  than  27.0.  The  dashed,  dotted,  and  solid  curves  indicate  the 
outflow  (eastward),  the  inflow  (westward),  and  the  net  transports. 


6)  at  150-m  depth  (around  25.0  crs  surface):  the  Ku- 
roshio  intrusion  (Fig.  10)  is  the  strongest  in  January, 
and  the  salty  tongue  (Fig.  6)  has  the  maximum  pene¬ 
tration  in  June.  Such  a  lag  (4-5  months)  is  approxi¬ 
mately  the  travel  time  for  the  Kuroshio  water  penetrating 
from  the  Luzon  Strait  into  the  southern  SCS.  For  a  rough 
estimation,  water  particles  in  the  winter  (December- 
February)  enter  the  Luzon  Strait  with  a  speed  around 
of  0.10  m  s_1  at  150-m  depth  and  travel  8.64  km  each 
day.  The  distance  between  the  Luzon  Strait  and  the 
southern  SCS  is  nearly  1200-1400  km.  It  will  take  ap¬ 
proximately  5  months  for  the  water  particles  traveling 
from  the  Luzon  Strait  to  reach  the  southern  SCS. 

8.  Conclusions 

The  goal  of  this  study  is  to  investigate  the  seasonal 
circulation  patterns  of  the  South  China  Sea  on  isopycnal 
surfaces  using  the  navy’s  unclassified  GDEM  data  on 
Yz°  X  Vi  resolution.  We  reformulate  the  P- vector  inverse 
method  for  determining  the  absolute  geostrophic  veloc¬ 
ity  on  the  isopycnal  surface.  The  circulation  pattern 
agrees  quite  well  with  earlier  observational  (Wyrtki 
1961;  Chu  et  al.  1997a,b;  Chu  et  al.  1998b)  and  mod¬ 
eling  studies  (Metzger  and  Hurlburt  1996;  Chu  et  al. 
1998c,  1999). 

(1)  The  GDEM  data  at  intermediate  depth  (150  m) 
shows  that  the  SCS  basin  is  dominated  by  a  cool  pool 
with  several  cold  centers  located  northwest  of  Luzon 
and  south  of  Dongsha,  near  the  South  Vietnamese  Bight, 


Xisha,  and  Liyue  Bank.  The  warm  water  appears  in  the 
vicinity  of  the  SCS  cool  water.  Two  other  warm  centers 
appear  southwest  of  Luzon,  and  northwest  of  Borneo. 
The  SCS  seasonal  thermal  variability  is  characterized 
by  strengthening/weakening  and  expanding/shrinking  of 
the  major  SCS  cool  pool.  The  basin-wide  cool  water 
expand  during  the  winter  and  shrinks  during  the  sum¬ 
mer.  In  December,  this  cool  pool  occupies  the  vast  area 
of  the  SCS  from  9°  to  20°N.  The  warm  tongue  penetrates 
northwestward  through  the  Luzon  Strait  from  the  Ku¬ 
roshio  and  then  stretches  southwestward  along  the  west¬ 
ern  boundary  of  the  SCS  with  a  very  narrow  width 
(100-200  km). 

(2)  The  GDEM  data  at  intermediate  depth  (150  m) 
shows  the  advance/retreat  of  the  salty  tongue  repre¬ 
sented  by  34.60  psu  isoline  from  the  Philippine  Sea 
through  the  Luzon  Strait.  This  salty  tongue  penetrates 
northwestward  right  after  passing  through  the  Luzon 
Strait  and  then  stretches  southwestward  along  the  west¬ 
ern  boundary  of  the  SCS.  The  salty  tongue  has  the  max¬ 
imum  (minimum)  penetration  into  the  SCS  in  June  (Oc¬ 
tober). 

(3)  The  subsurface  (crg  =  23.0  kg  m~3)  circulation 
has  a  strong  seasonal  variation  with  a  basin-wide  cy¬ 
clonic  gyre  in  the  winter  and  a  weak  anticyclonic  gyre 
in  the  summer.  Warm  and  salty  Kuroshio  water  enters 
the  SCS  through  the  Luzon  Strait  evidently  in  October- 
March,  the  transition  to  and  during  the  winter  monsoon 
season.  A  cyclonic  mesoscale  eddy  splits  from  the  gyre 
near  Dongsha  in  February,  strengthens  afterwards,  and 


September  2000 


CHU  AND  LI 


2437 


becomes  evident  through  the  whole  summer  monsoon 
season.  The  Kuroshio  intruded  flow  is  the  north  or  the 
northeast  fringe  of  the  DS  eddy.  The  maximum  speed 
reaches  0.2  m  s  1  in  April.  The  Dongsha  cyclonic  eddy 
is  undetectable  from  November  to  January.  During  the 
summer,  the  northward  western  boundary  current  occurs 
from  the  south  to  the  Vietnamese  coast  with  a  much 
weaker  speed  (less  than  0.1  m  s_1)  in  the  summer  than 
its  counterpart  (southward  western  boundary  current)  in 
the  winter.  The  western  boundary  current  splits  into  two 
currents  at  10°N.  The  offshore  current  recirculates  and 
forms  a  mesoscale  anticyclonic  eddy  near  Natuna  Island. 
The  major  feature  of  the  circulation  at  the  intermediate 
level  (cre  =  26.2)  is  the  persistence  of  the  basin-scale 
cyclonic  gyre.  There  is  no  seasonal  reversal  of  the  west¬ 
ern  boundary  current.  The  gyre  is  evident  and  strong  in 
the  winter  and  weakens  in  the  summer.  As  spring  ap¬ 
proaches  (April-May),  the  western  boundary  currents 
(southward)  near  the  S  VB  weakens.  This  “broken”  gyre 
feature  continues  in  the  whole  summer  (June-August). 
The  western  boundary  currents  (southward)  intensifies 
in  the  fall. 

(4)  We  computed  monthly  total  and  layered  (between 
two  cre  levels)  volume  transports  through  the  Luzon 
Strait  along  120°45'E  longitude.  The  outflow  is  much 
weaker  than  the  inflow  for  the  total  transport  and  the 
most  layered  transports  except  the  deep  layer  (crg  > 
27.025),  where  the  outflow  is  dominant.  The  total  trans¬ 
port  is  negative  (inflow)  all  year  round  with  a  minimum 
value  of  — 13.7  Sv  in  February  (strongest  intrusion)  and 
a  maximum  value  of  —1.4  Sv  in  September  (weakest 
intrusion)  with  the  annual  transport  of  —6.5  Sv.  These 
values  may  change  when  different  cross  section  is  used. 

(5)  Current  ocean  numerical  models  are  initialized 
from  either  observed  or  climatological  temperature  and 
salinity,  as  well  as  zero  velocity.  This  leads  to  an  im¬ 
balance  between  the  velocity  and  density  fields  during 
the  spinup  period.  If  we  use  the  geostrophic  velocity, 
as  the  initial  condition,  the  velocity  and  density  fields 
are  balanced.  Thus,  the  inverse  method  proposed  here 
provides  a  better  initialization  for  the  isopycnal  models. 

(6)  Two  important  parameters,  volume  transport 
through  the  Taiwan  Strait  and  wind  stress,  are  not  in¬ 
cluded  in  this  study.  Exclusion  of  the  volume  transport 
through  the  Taiwan  Strait  is  due  to  the  difficulty  of  the 
inversion  on  the  isopycnal  surface  in  shallow  depths. 
Without  the  wind  stress,  the  circulation  presented  here 
is  only  the  geostrophic  flow.  Chu  et  al.  (1998c)  showed 
that  the  wind  stress  largely  affects  the  seasonal  vari¬ 
abilities  (especially  in  the  circulation)  in  upper  layer  (30 
m)  of  SCS.  Thus,  if  the  wind  stress  is  considered,  the 
seasonal  variabilities  at  the  sub-surface  (cr8  =  23.0  kg 
m  !)  are  expected  to  be  larger.  Inclusion  of  wind  stress 
in  the  velocity  inversion  is  one  of  our  on-going  projects. 
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